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The use of the absolute measure of local chirality is championed since it has a uniform distribution 
for randomly reshuffled chiral components so that any deviations from uniformity in the associ- 
^ — , | ated "X-distribution" are directly attributable to QCD-induced dynamics. We observe a transition 

in the qualitative behavior of this absolute X-distribution of low-lying eigenmodes which, we 



propose, defines a chiral polarization scale of the QCD vacuum. 
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1. Introduction 

Important properties of the QCD vacuum can be studied indirectly by looking at the character- 
istics of individual low-lying overlap fermion eigenmodes, which provide a natural way to filter out 
UV fluctuations. The space-time chiral structure of the eigenmodes reflects the space-time topolog- 
ical structure of the underlying gauge fields. Properties of the "X-distribution"[[I]], the probability 
distribution of the local chiral-orientation parameter, stirred a debate [||] whether this offered evi- 
dence in favor of, or against, models of the QCD vacuum. We revisited [^, ||] the use of the local 
chiral-orientation parameter and found that most of the qualitative features of the distribution were 
a kinematical effect which was removed with the use of an improved and absolute measure of lo- 
cal chirality. Here we flesh out the construction of the "absolute X-distribution" and proceed to 
analyze the striking change in its behavior as the energy scale of the eigenmodes is scanned [0]. 

2. Local Chirality and the X -Distribution 

The local-chirality parameter measures the tendency of a low-lying Dirac eigenmode, \j/ = 
Wl + Wr to be left or right handed. Sampling the space-time values of a particular set of eigen- 
modes, say at a particular energy scale, for a set of configurations yields (a statistical approximation 
of) a base probability distribution ^ (<7i , #2) where qi = \yl\ and q2 = \ Yr\. A scatter plot of such 
a distribution for a particular set of modes is shown in the left panel of Figure [I]. We can elucidate 
the tendency for polarization by integrating over the radial coordinate, p, and investigating the de- 
pendence on the polar angle, (p. It is convenient to symmetrize the angular variable <p € [0,7t/2] 




Figure 1: Left: Dynamics governing the polarization properties of two lowest modes for ensemble E\ 
(defined in Table 1). Right: The associated "uncorrelated" dynamics; solid gray lines separate the quadrant 
into 10 sectors, each containing 10% of the population describing 

by defining the "reference polarization coordinate" x = ^(p — 1 so that x € [— 1,+1] and such that 
a value of x = +1 (x = — 1) corresponds to purely right handed (left handed) bispinors [ffj. Thus, 
in the literature [^], histograms which have a "double-peak," that is two peaks near x = ±1 in a 
distribution symmetric about x = 0, have been interpreted (prematurely as we will see) as revealing 
that the eigenmodes display chiral polarization. The amount of peaking can be enhanced with a 
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different choice of "polarization function," that is, using functions of x, rather than x itself. These 
polarization functions can be classified as belonging to one of several families, such as 



X c (x;a) =sgn(x) \x\ a ; X R (t;a) 



1 



t a + \ 



X G (t;a) = - tan 



Ht a ) 



1 



(2.1) 



where a > and t = tan(<p). X c (x; 1) = x = X G (t; 1) correspond to the original chiral orientation 
parameter X G (t;2) and X R (t;2) are alternatives that appeared in the literature [0|. All of 
these choices of polarization functions share the same crucial features of the simplest choice of the 
reference polarization function X(x) =x, namely they are odd, monotonically-increasing functions. 





• a=4 







0.0 
X 
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Figure 2: Possible ^-distributions associated with fixed dynamics can have qualitatively different behavior. 
Here we select two polarization functions from family X c (left), and two from family X R (right). In this and 
subsequent figures, the data have error bars but they may to be too small to be resolved. 
The marginal distribution 



P(X) = J dq x J dq 2 & b {q x ,q 2 )d(x-X{q 1 ,q 2 ) 



(2.2) 



where X is a generic independent variable parametrizing the range of polarization functions, i.e. 
X 6 [—1,+1], is called the "X -distribution" [jl]|. It depends not only on the dynamical base proba- 
bility S^b^qxA'i) but also on the polarization function X(x) used to measure it. As Fig. ||| indicates, 
the qualitative features of the X-distributions depend very strongly on the choice of the polariza- 
tion function, and one might be led into declaring, perhaps falsely, that the distribution is highly 
polarized or highly unpolarized depending on ones arbitrary choice. 

Indeed, most of the qualitative appearance of an X-distribution is a result of "kinematics" (or 
"phase space") and is not due to QCD correlations at all [ph. To demonstrate this, we compute 
the X-distribution arising from i/^(y) and i/^(y), where y is the spacetime coordinate, and then 
randomly reshuffle the fields YLran(y) = Yi{yran), where y ran is a random permutation, using an in- 
dependent random reshuffle for Any QCD-dynamically-induced correlation between Vz,(y) 
and y/#(y) at site y is destroyed by this procedure. Nevertheless, the correlated and uncorrected 
scatter plots of the base probability distributions are very similar as is seen in Fig. |l[ Accordingly, 
when we re-compute the X-distribution, the randomized uncorrected X-distribution looks very 
similar to the original as seen in Fig. ||[ 



3. The Absolute X-Distribution 

To expose the true correlation we need to "subtract" the phase space background, or more 
precisely measure the (QCD-dynamically-generated) correlated distribution relative to the (ran- 
domized) uncorrected one. 
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Reference X -distribution for ensemble E.\{24 A ) 



Reference X— distribution for ensemble £\|(24" 1 ). uncorrected 
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Figure 3: The references-distribution for "lowest" and "highest" modes in ensemble £4 (defined in Table 1) 
is shown in the left panel. The associated distribution with statistically independent left-right components 
("uncorrelated") is shown in the right panel. 

The following procedure makes a differential comparison of polarization in the dynamics 
&b(qi,<l2) to that of uncorrelated probability distribution ^^(91,92) = P^Opfe)' where p{q\) = 
J -i dqi&bifHtqi) 1S a marginal distribution: A ray that passes through the origin is specified by 
its reference polarization coordinate x, as shown in Fig. [I]. Determine the fraction of population 
contained between the gi-axis and this ray, separately for &b{qi,qi) an d for ^^(qi, qi)- These 
fractions are the cumulative probability functions S r (x) = J^ l dx'P(x') and S"(x) = J^dxfP^x'). 
Eliminating x, every ray is represented by a point in the cumulative distribution plane (S^,S r ) with 
the set of all such points forming a curve starting at (0,0) and ending at (1,1). If the the two 
dynamics in question have identical polarizations the graph will be a straight line. On the left side 
of Fig. || we plot the result of this construction for data displayed in Fig. [j] The correlated and 
uncorrelated dynamics have very similar polarizations since the graph is close to being linear. To 
obtain a differential comparison, and to better see the differences, we compute the slope of the 
cumulative polarization graph and show it on the right side of Fig. [| As one can clearly see now, 
the correlated dynamics exhibits a small excess of polarization with respect to the uncorrelated one 
near the extremal values. This is a plot of the "absolute X -distribution" 



Pa 



1 dS r 

YdSi 



(3.1) 



except for the additional rescaling S" — > X = 2S" — 1 required since X € [— 1, + 1] while the un- 
correlated cumulative probability S" 6 [0,+l]. This construction is manifestly independent of the 
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Figure 4: Steps in the construction of absolute X-distribution for data shown in Fig. 1. See the discussion 
in the text. 
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choice for polarization coordinate, with different possibilities corresponding to different parametric 
representations of the same curve in the (S",S r ) plane. 

Fig. H shows the absolute X-distribution Pa(X) for the lowest two non-zero pairs (left panel) 
and the highest two (of those measured) from ensemble £5 (defined in Table |T]). In comparison is 
shown the uncorrected X-distribution, P"(x), constructed using the reference polarization coordi- 
nate x. The shape of Pa{X) is determined by the correlation induced by QCD dynamics, while the 
shape of P"(x) is determined completely by the kinematics. From the pair, Pa(X) and P"(x), one 
could reconstruct the correlated (and uncorrected) X-distributions for any choice of polarization 
function X(x), such as those chosen in [gj. 



Lowest 2 modes distributions for ensemble £ 5 (16 4 ) 

— • — absolute 
— ■ — uncorrected 



Highest 2 modes distributions for ensemble E 5 (16 4 ) 



absolute 
uncorrelated 



Figure 5: Absolute (dynamical) and uncorrected (kinematical) X-distributions for gauge ensemble £5 
(defined in Table 1 ) using the reference polarization coordinate for the latter. Results for lowest two non-zero 
pairs are shown on the left while the highest two pairs are shown on the right. 



4. Convexity or Concavity of Absolute X-Distributions 

Now we look at the X-distributions for the lowest and highest modes (of those measured low- 
lying modes) for our finest (24 4 ) ensemble (£4). We have a total of four ensembles on lattices at 
various lattice spacings, with fixed volume, and a fifth to check for finite-size effects. These are 
listed in Table |l| 



Ensemble 


Size 


config 


Volume 


Lattice Spacing 


A MAX 
yv LOW 


aAVE 
iv LOW 


aMIN 
HIGH 


aAVE 
HIGH 


Ei 


8 4 


100 


(1.32fm) 4 


0.165fm 


449 


226 


1956 


1980 


E 2 


12 4 


97 


(1.32fm) 4 


O.llOfm 


407 


169 


1711 


1735 


£3 


16 4 


99 


(1.32fm) 4 


0.0825 fm 


304 


142 


1513 


1553 


£4 


24 4 


96 


(1.32fm) 4 


0.055 fm 


344 


136 


1338 


1366 


E 5 


16 4 


99 


(1.76fm) 4 


O.llOfm 


162 


58 


1087 


1123 



Table 1: The summary of five ensembles used in overlap eigenmode calculations. The right side of the table 
describes some properties of the spectra (in MeV) with A^qw denoting the average magnitude of lowest 
near-zero eigenvalue over the ensemble and A^)q H denoting the same for highest eigenvalue. A^ovy is the 
magnitude of the maximal lowest eigenvalue, and AjjjQ H the magnitude of the minimal highest eigenvalue. 



Figure |3| (left panel) shows the (reference) X-distribution for the lowest pair of eigenmodes and 
for the highest pair (of the 55 measured). "Subtracting" from this the corresponding X-distribution 
for randomized data leaves the absolute X-distribution, Pa(X), shown Fig. |[ We see that much 
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of the peaking in the reference X-distribution has been removed, by the comparison to the uncor- 
related version. This indicates that most of the peaking was kinematical rather than dynamical in 
the original reference distribution and that to draw any conclusions from this distribution would be 
misleading at best. On the other hand, for the absolute X-distribution, any deviations from unifor- 
mity are attributable to correlations induced by QCD dynamics. There is a relatively small but clear 
effect exhibited by the lowest modes which have a convex shape and the higher modes which have 
a concave X-distribution. That is to say, the lowest modes exhibit an enhanced chiral polarization 
while the higher modes exhibit a suppression, which we regard as a novel and tantalizing discovery. 



Absolute X-distribution for ensemble £ 4 (24 4 ) Absolute X-distribution for ensemble £ , 4 (24 4 ) 
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Figure 6: A plot of the absolute X-distribution for the "lowest" (left panel) and "highest" (right panel) 
modes for our finest ensemble. 

5. Chiral Polarization Scale 

From Figs. [6| and [7] we see that the shape of the absolute X-distribution curve is qualita- 
tively different for the lowest several modes (convex) versus higher modes (concave). We seek 
the scale A^ at which the transition from convex to concave occurs. To do this we first define, as 
a figure of merit, a dynamical correlation coefficient in terms of the first moment of the absolute 
X-distribution. 

C A (X) = -l + 2 J dX\X\P A {X) 

for which Ca > for enhanced polarization, Ca < for suppressed polarization, and Ca = for 
neutrality as for the case of statistical independence. 

Fig. H (left panel) shows that the dynamical correlation coefficient decreases monotonically 
from positive values at low modes to negative values at higher modes. A linear interpolation de- 
termines the scale At, at which it crosses zero, at which point there is no dynamical tendency for 
enhancement or suppression of local chirality in comparison to the associated X-distribution for 
randomized, uncorrected left-right components. 

We repeated the determination of this scale for each of four lattice spacings in order to look 
for a continuum limit. Fig. || (right panel) suggests that Aj(V,a) has a finite continuum limit at 
fixed volume. A fit to the form c\ +C2t? 4 , with the coarsest lattice excluded from the fit, is drawn 
to guide the eye. Further studies must address whether this "chiral polarization scale" survives 
in the infinite-volume limit, but our check of the finite volume effect at a single lattice spacing is 
encouraging. 
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A =350 MeV 



A=420 MeV 
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Figure 7: Absolute ^-distributions for gauge ensemble £4 with changing scale A of the modes. Note the 
transition from convex to concave (positive to negative correlation) between 470 MeV and 530 MeV. 

Continuum limit of the transition poinl 



Dynamical correlation coefficient for ensemble E 4 




Figure 8: Left panel: Determination of the chiral polarization scale, Ar. for the finest lattice. Right panel: 
Dependence of the chiral polarization scale, At(V,o) on the lattice spacing at finite volume. 

6. Summary 

X-distributions are detailed probes of QCD dynamics. The use of previous definitions in the 
literature is misleading, since the distributions are dominated by kinematical effects. It is imperative 
to use the absolute X-distribution as it removes kinematical effects leaving correlations induced by 
QCD dynamics. We have discovered that the absolute X-distributions of the lowest several Dirac 
eigenmodes are convex (polarized) and those of higher modes are concave (anti-polarized). We 
have determined the scale Aj of the transition; its continuum limit is finite (at least for finite 
volume) and, it is proposed, can define a "chiral polarization" scale of QCD. Complete details are 
in reference 101. The absolute X-distribution was first discussed in [Si. 
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